Decay process accelerated by tunneling in its very early stage 
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We examine a fast decay process that arises in the transition period between the Gaussian and 
exponential decay processes in quantum decay systems. It is usually expected that the decay is 
decelerated by a confinement potential barrier. However, we find a case where the decay in the 
transition period is accelerated by tunneling through a confinement potential barrier. We show that 
the acceleration gives rise to an appreciable effect on the time evolution of the nonescape probability 
of the decay system. 
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Quantum decay processes such as the nuclear alpha 
decay are usually well described by means of the expo- 
nential decay law; see, for example, [|l|, ^, ^. Theoret- 
ically, however, deviations from the exponential law are 
expected in the beginning and also toward the end of 
a decay process It is understood that the quantum 
decay process in general proceeds through three differ- 
ent stages; initial, intermediate and final. The initial 
stage is characterized by the Gaussian law, the interme- 
diate stage by the exponential law, and the final stage by 
the power law. The decay speeds in the initial and final 
stages are smaller than that in the intermediate stage. 
In particular, the slow decay process in the Gaussian pe- 
riod leads to the possibility of the quantum Zeno effect 
Q , in which the decay process is decelerated by repeated 
measurements. 

Many years ago, in his one-dimensional model analy- 
sis of a decaying system. Winter found that the speed of 
the decay process exhibits an irregular behavior in the 
transition period between the Gaussian and exponential 
periods Very recently Dicus et al. reexamined the 
same system, which consists of a particle which is initially 
confined in a region and leaks out by tunneling through 
a delta- function potential barrier ^ . In the irregular de- 
cay process, the decay speed can be larger than that in 
the exponential region. Such fast decay in the transition 
period is interesting in the sense that it may give rise to 
acceleration of the decay process by repeated measure- 
ments, which is the so-called anti- (inverse-) Zeno effect 

The purpose of this report is to investigate the details 
of the fast decay process in the transition period between 
the Gaussian and exponential periods |9j. In particular, 
we explore how the decay process in the transition period 
depends on the strength of the confinement potential bar- 
rier with a finite potential width. It is usually expected 
that the decay is suppressed as the confinement becomes 
stronger because the potential barrier suppresses the time 
evolution of the wave function. However, we will illus- 
trate a case where the decay in the transition period is 
accelerated by tunneling through a confinement poten- 
tial barrier. Further, we show that the acceleration gives 



rise to an appreciable effect on the time evolution of the 
nonescape probability of the decay system. 

We consider a model in one dimension with a potential 
that consists of an infinite repulsive wall at x = and a 
repulsive square barrier atl<2:<l + w, 

( for < a; < 1, 
V[x) = I h iov Kx <l + w, (1) 
I for 1 + w < X, 

where h and w are the height and width of the potential, 
respectively. 

We assume that a particle of mass m is initially con- 
fined within the potential barrier and it leaks out in time. 
The time-evolution of the system is determined by the 
time-dependent Schrodinger equation 

dMx,t) , , 

'^^ = H4>{x,t\ H^-— + V{x) (2) 

where the units are such that h — I and 2m 1. 
We solve Eq. (H) numerically using the implicit solu- 
tion method jlO|, with the unitary time-evolution opera- 
tor f7 = (1 - iH6t/2)/{l -f iHSt/2), where St is the time 
mesh. For the range of x, we take [0,500]. In order to 
suppress reflection of the wave function at the boundary 
{x = 500) of the model space, we assume an imaginary 
potential at x > 490. 

For the initial wave function of the particle which is 
confined in the region of < a; < 1, we assume 

0) = 72 siuTra;. (3) 

It is understood that 0(a;,O) = for x > 1. The 
wave function leaks out in time by tunneling through the 
potential barrier V{x). The energy expectation value 
< > at i = is TT^. In numerical illustrations 
we choose the height of the potential barrier such that 
h > n^. In this sense we regard the process as a tunneling 
decay process. 

We define probability P{a,t) for the particle being in 
the interval < x < a aX t hy 

P{a,t)= f dx\^ix,t)\^. (4) 
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FIG. 1: The time-evolution of the nonescape probabihty 
P(4, t) for w = 0.6 and h = 10. The units are such that h = 1 
and 2m = 1. 



FIG. 2: The time-evolution of ^(4, t) for a fixed width w = 
0.6. The solid, dashed and dot-dashed lines show g(4, t) for 
the heights h — 10, 20 and 30, respectively. The units are 
such that h = 1 and 2m = 1. 



If we set a = 1 + w, P{l + w,t) represents the probability 
that the particle is confined by the potential barrier. For 
w we consider a few different values. In order to compare 
the results for different values of w, we set a such that 
a > 1 + Wmax- Throughout this paper, we take a = 4. 
The P{a,t) is a kind of the nonescape probability, which 
is the probability that the particle has not escaped from 
the potential by time t 
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Next we introduce function g{a,t) defined by 



g{a,t) 



dP{a,t)/dt 



(5) 



If probability P{a,t) decays exponentially, that is, 
P{a,t) cx e~^*, then g{a,t) is independent of time. 



if 



g{a,t) 



'1- 



If P{a,t) obeys the Gaussian decay law, that is, if P cx 
e~* Z"^, g{a,t) is proportional to t: 



g{a,t) = -2t/T. 



(6) 



Thus, from the t dependence of 17(0, t) we can see how 
well P{a, t) obeys the exponential law or the Gaussian 
law. 

Figure 1 shows the nonescape probability P{a,t), 
where the potential height and width are taken as /i = 10 
and w — 0.6, respectively. In the initial period < t < 
0.3, the decay is extremely slow. The period should cor- 
respond to the Gaussian decay process. In the period of 
t >2, the system is subject to the exponential decay law. 
Between the two stages, there is a period (0.3 ^ t < 2) 
in which the decay process is neither Gaussian nor ex- 
ponential. This is the transition period. In the period 
the decay of the nonescape probability seems to be much 
faster than that in the exponential period. 

The function g{4:,t) calculated from P(4,t) of Fig. 1 
is shown by the solid line in Fig. 2. In the period of 
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FIG. 3: The time-evolution of (?(4,t) with a fixed height 
h — 15. The solid, dashed and dot-dashed lines are the t;(4, t) 
for w = 0.6, 0.8 and 1.8, respectively. The units are such that 
h — 1 and 2m = 1. 



t > 2, g{4:,t) is almost constant, which means that the 
decay process is exponential. In the period of < i < 0.3, 
5(4, t) is not exactly proportional to t, i.e., the decay pro- 
cess in the initial stage slightly deviates from the Gaus- 
sian decay law. However, the decay speed in this period 
is still smaller than that in the exponential period. Cor- 
responding to the rapid decay process seen in P{t) of 
Fig. 1, the maximum of |(7(4,t)| is obtained in the tran- 
sition period {t ^ 0.6). The g{4:,t) starts from zero at 
t = 0. Thus, the quantum Zeno effect is possible when 
we repeat measurements with a sufficiently small time 
interval. On the other hand, the anti-Zeno effect is pos- 
sible by repeated measurements only when the net decay 
rate of the fast decay process in the transition period is 
large compared to that of the slow decay process in the 
Gaussian period (see footnote p3{). 
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FIG. 4: The time-evolution of (;(4,t) for a fixed width w — 
0.2. The dashed, dot-dashed and dotted lines are the (?(4, t) 
for the widths h = 10, 20 and 30, respectively. The solid line 
shows the g{4:,t) for no potential barrier. The units are such 
that h — 1 and 2m — 1. 



In order to see how the decay process in the transi- 
tion period depends on the strength of the potential bar- 
rier, we examine (7(4, t) for various potential heights and 
widths. The dashed and dotted-dashed lines in Fig. 2 
show the g(4, t) for h — 20 and h = 30, respectively with 
a fixed width w — 0.6. The decay process in the Gaus- 
sian period < t < 0.3 does not depend strongly on 
the potential height. This is because in the initial stage 
the higher energy components of the initial wave func- 
tion relative to the potential height contribute mainly to 
the decay of the system. On the other hand, the decay 
speed in the exponential period becomes much smaller 
as the confinement becomes stronger. The fast decay in 
the transition period depends strongly on the potential 
height. It tends to be suppressed as the confinement be- 
comes stronger. 

In Fig. 3, we show ^(4, t) for various potential widths 
and a fixed height h = 15. The decay process in the Gaus- 
sian period < t < 0.3 does not seem to depend strongly 
on the potential width. On the other hand, the decay 
rate of the exponential period becomes much smaller as 
the potential width becomes broader. This is due to the 
increase of the confinement strength. The decay speed 
in the transition period becomes larger as the potential 
width becomes narrower. 

As we have shown, the speed of the fast decay pro- 
cess becomes smaller as the potential barrier becomes 
stronger. Thus, we might guess that the fastest decay 
will be obtained in the decay process with no potential 
barrier. However, as we will show in the following, po- 
tential barriers with appropriate widths and heights can 
accelerate the decay process. 

In Fig. 4, the solid line shows g{4:,t) calculated with 
no potential barrier. The dashed, dot-dashed and dotted 
lines exhibit 5(4, t) calculated with the potential barriers 
with h = 10, 20 and 30 and fixed width w = 0.2, respec- 
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FIG. 5: The time evolution of the nonescape probabilities 
P(4, t). The solid and dashed lines are the P(4, for no 
potential barrier and for the potential barrier with w = 0.2 
and h = 10, respectively. The units are such that h = 1 and 
2m = 1. 



lively. It should be noted that w = 0.2 is much thinner 
than those used in Figs. 2 and 3, and therefore, in this 
case, the confinement is very weak compared with that 
in Figs. 2 and 3. In the Gaussian period, the potential 
dependence of the decay speed is not very appreciable. 
However, in the transition region, the decay speed be- 
comes faster as the potential height becomes lower. For 
h = 10, the maximum decay speed at i ~ 1.8 exceeds 
that for no potential barrier at t ^ 0.7. 

As shown in Fig. 5, such an acceleration of the decay 
speed by tunneling gives rise to an appreciable difference 
in the time evolution of P{A,t). The solid and dashed 
lines are the P(4, t) for no potential barrier and for the 
potential barrier with w = 0.2 and h = 10, respectively. 
The nonescape probability for /i = 10 becomes smaller 
than that for no potential barrier at t ~ 1.5. At this 
time the residual nonescape probability is still about ten 
percent. In this sense, the effect of this acceleration can- 
not be ignored. On the other hand, for h = 20, at the 
time region in which the nonescape probability becomes 
smaller than that for no potential barrier, the residual 
nonescape probability is negligibly small. 

The fluctuations in the behaviors of g(4, t) indicate 
that the decay processes are still in the transition pe- 
riod from the Gaussian to exponential period. However, 
as shown in Fig. 5, the decay process has been almost 
completed before t = 4. Therefore, even if the decay pro- 
cess proceeds to the exponential period eventually, the 
exponential decay has no importance in this case. The 
stability of quantum system is usually characterized by 
the magnitude of the imaginary part of the pole that 
gives the inverse of the lifetime in the exponential period 
IQ, 1^, ^ . The result that we have shown implies that such 
a pole analysis may not be effective for highly unstable 
systems Q. 

In Fig. 4, we have seen that the acceleration by tun- 
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FIG. 6: The time-evolution of (/(4, i) with a fixed height 
h = fO. The dashed, dot-dashed and dotted lines are the 
g{4:, t) for the heights w — 0.2, 0.4 and 0.6, respectively. The 
solid line shows the 3(4, t) for no potential barrier. The units 
are such that h = 1 and 2m — 1. 

neling can be obtained for h < 10. Next, we investigate 
the acceleration with a fixed height. We examine the 
time evolution of g{4:,t) for ft, = 10 |l^. The dashed, 
dot-dashed and dotted lines in Fig. 6 exhibit g{4, t) cal- 
culated for the potential barriers with w = 0.2, 0.4 and 
0.6 and fixed height h = 10, respectively. One can see 
that the decay rate becomes larger for thinner potential 
widths. For w = 0.2, the maximum decay speed ex- 
ceeds that for no potential barrier. Thus, one sees that 
the acceleration of the decay speed by tunneling can be 
obtained when the strength of the confinement by the 



potential barrier is sufficiently weak. In our illustrations, 
for h <10 and w < 0.2, the accelerations are remarkable. 

Finally, we mention that the g(4, t) represented by dot- 
ted line in Fig. 4 {w — 0.2 and h = 10) takes positive 
values around t ~ 3.5. Recall that the probability cur- 
rent J (a, t) is related to the nonescape probability P(a, t) 
by j{a,t) — —dP{a,t)/dt ~ —g{a,t)P{a,t). This means 
that, if g{a,t) is positive, j{a,t) is negative. However, in 
|l6j the negative currents were obtained at very late 
time region after the exponential decay period. Our re- 
sult implies that in a highly unstable quantum system 
the negative current can occur even at the initial stage 
after the Gaussian period. 

We have investigated the fast decay process in the tran- 
sition period between the Gaussian and exponential de- 
cay processes. In most cases of the tunneling process, 
the decay speed becomes smaller as the potential barrier 
becomes stronger. As a special case, we have found that 
the fast decay process can be remarkably accelerated by 
tunneling through potential barriers with appropriately 
small widths and heights. A detailed analysis of the ac- 
celeration of the fast decay process by tunneling is a fu- 
ture project. 
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